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Abstract 

The neutrino long wavelength (just-so) oscillation is reconsidered as a solution to 
the solar neutrino problem. In the light of the presently updated results of the four 
solar neutrino experiments, the data fit in the just-so scenario substantially improves 
and becomes almost as good as in the MSW scenario. Surprising result of our analysis is 
that best fit is achieved when the oscillation occurs only between two neutrino states: 
switching on the oscillation into third neutrino increases the x 2 value. Namely, we 
consider the vacuum oscillation scenario in the three-neutrino system (4 parameters) 
and find out that the x 2 minimum is always achieved in the two parameter subspace 
in which actually only two neutrino states oscillate. This holds in the framework of 
any solar model with relaxed prediction of the various neutrino fluxes. The possible 
theoretical implications of this observation are also discussed. 
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There are strong arguments to believe that the Solar Neutrino Problem (SNP), i.e. 
the deficit of the solar neutrino fluxes indicated by four operating experiments []T|, ^ 
U as compared to the predictions of the Standard Solar Models (SSM) || ||, cannot be 
explained by the nuclear /astrophysical reasons. Namely, the results of different detectors 
cannot be reconciled among each other by varying the SSM parameters as are the central solar 
temperature, nuclear cross sections etc. 0. The problem is essentially related to the peculiar 
energy dependence needed for the suppression of different solar neutrino components, so that 
the intermediate energy 7 Be neutrinos are to be killed more than the lower energy (pp) or 
higher energy ( S B) neutrinos. This points that the SNP is rather due to the "non-standard" 
neutrino properties. The most natural and plausible possibility is the oscillation of the solar 
u e into another type of neutrino v x . 

The neutrino oscillation picture can provide the necessary energy dependence in two 
regimes, which are known as the Mikheyev-Smirnov-Wolfenstein (MSW) || and the long 
wavelength vacuum oscillation (so called just-so) ][| scenarios. The MSW resonant conversion 
in matter offers a natural possibility of selective strong reduction of the 7 Be neutrinos and 
thus appears to be the most attractive and elegant solution to the SNP. It provides a very 
good fit of the experimental data, implying the mass range Sm 2 ~ 10~ 5 eV 2 and small mixing 
angle, sin 2 29 ~ 10~ 2 [|1^]. The just-so scenario, with the oscillation wavelength comparable 
to the sun-earth distance, needs Sm 2 of about 10~ 10 eV and large mixing angle, sin 2 2(9 ~ 1 
13|, |T4f, which parameter range can be naturally generated by non-perturbative quantum 



gravitational effects |I| [TT 



Both the MSW and just-so mechanisms are theoretically motivated^] and so far the ex- 
perimental data do not allow to discriminate them. Hence, both these scenarios should 
be considered as realistic candidates to the SNP solution. However, the peculiar "just-so" 
predictions as are the seasonal time variation of the neutrino signals and the specific defor- 
mation of the original solar neutrino spectra 0, 14], will allow to test the long wavelength 
oscillation scenario at the future real time detectors like Superkamiokande, SNO, Borexino, 
etc., and to discriminate it from the MSW picture. 

In the present paper we present an updated analysis of the just-so scenario. The earlier 



analysis [T3L 14] has shown that the experimental data fit in this scheme is acceptable, but 



somewhat worse than that of MSW picture. However, after publication of these papers the 
experimental data have been changed. In particular, the result of the Homestake experiment 
has been recalibrated [HJ and the new data of GALLEX experiment became available ||. 
On the other hand, the SSM predictions have been reconsidered by taking into account 
the helium and heavy element diffusion ||. The changes are small and rather confirm the 
stability of both the experimental and theoretical results. Nevertheless, as we show below, 
these small changes make the "just-so" data fit substantially better so that it becomes almost 
as good as in the MSW picture. 

Another issue which we address in this paper is the following. Up to now, in the literature 
all the analysis of the vacuum oscillation scenario was performed within the simplified case 
of two neutrino (2z/) system, which case employs only two parameters: the v e — v x mixing 



1 In the language of the Standard Model (SM) the neutrino masses emerge through the higher order 
operators of the type -jj{liClj)HH , where li (i — 1,2,3) and H are respectively the lepton and Higgs 
doublets and M is some regulator scale. (For example, the famous "seesaw" scenario effectively reduces to 
these operators with M of the order or the right-handed neutrino mass). Then the neutrino mass range 
needed for the just-so scenario corresponds to the Planck scale, M ~ 10 19 GeV, whereas the MSW scenario 
requires M to be of the order of the supersymmetric grand unification scale, M ~ 10 16 GeV. 



angle 9 and the mass square difference Sm 2 between two neutrino eigenstates. However, the 
nature knows that there exist three neutrino species and if their masses are originated by the 
same mechanism, then the relevant case should be the oscillation picture in the 3v system, 
in which case the solar neutrino data are determined by four parameters. Here explore the 
full parameter space for the 3 neutrino just-so oscillations. As we will see, surprisingly the 
introduction of additional parameters dos not improve the data fit and the best description 
of the existing data is essentially provided by the 2v oscillation, i.e. by the case in which the 
oscillation into the third neutrino is decoupled. We briefly discuss the implications of this 
observation for the structure of the neutrino mass matrix. 



We perform our analysis in the spirit adopted earlier for the MSW picture in refs. [[L2 



and for the just-so scenario in our previous paper []l4j . Namely, instead of studying the 
just-so scenario in the context of one particular SSM with including in the x 2 analysis the 
theoretical uncertainties, we prefer to scan various realistic SSMs. We take as a reference 
model the recent Bahcall and Pinsonneault (BP) SSM || without the underlying theoret- 
ical uncertainties, while the variations of the SSM parameters are taken into account by 
parametrizing various neutrino fluxes as <p k = fk<fit =8 B, 7 Be, 13 N, 15 0, pep and pp), 
where the subscript '0' denotes the values predicted by the BP SSM || and the factors /& 
account for theoretical uncertainties. Thus, the cases of other SSMs || will be effectively 
recovered by relaxing <p k . In particular, for a SSM characterized by a given set of the 
signals of the radiochemical detectors are expected to be 

R s c f M = 7.36/b + 1.24/se + 0.22/ pep + 0.11/W + 0.37 f (SNU) 

R s G s a M = 69.7 f pp + 16.1/b + 37.7 f Be + 3.0f pep + 3.8f N + 6.3/ (SNU) (1) 

as compared to the experimental rates R^f = 2.55±0.25 SNU [|TJ and Rq% = 74±8 SNU (for 
the gallium data we use the weighted average of the GALLEX result Rgallex = 77 ± 8 ± 6 
SNU § and the SAGE one R S age = 69 ±11 ±6 SNU gj). Thus, the observed signal in each 
detector can be given as a ratio of the measured rate to that is expected in a given SSM, 
Z a = Rl xp /R^ SM (a = CI, Ga), and it is a function of the factors On the other hand, 
the Kamiokande result |2j implies that the ratio of the experimental rate to the prediction 
of SSM with given fs is 



R™ v 1 



— (0.45 ±0.08) (2) 



TDjO-M f 

K K JB 

where obviously the case /b = 1 corresponds to the BP SSM 0. 

Several numerical studies have shown that in practice the effects of independent varia- 
tions of the metal fraction, the opacities, the astrophysical factor S pp and the age of the sun 
are well reproduced by variations of the central solar temperature T c according to the power 
lows [r5|: 



(3) 



Thus, we consider /# = (T C /T C °)^ B as a free parameter while for other coefficients we cor- 
respondingly obtain /& = /| fc , where = (3k/ (3b- Following ref. |^5| , we take £ p 



-0.034 ±0.006, i Be = 0.47 ±0.08, 6v = 1.0 ±0.4, £ = 1.2 ±0.4 and £ pep = ±0.1.0 In this 

2 In the following computations we take the central values neglecting the theoretical uncertainties. These 
uncertainties are indeed small for the pp and Be neutrinos. For the CNO components these are bigger but 
the contribution of the latter to the signals (GT) are not substantial. 



way, the case of other SSM's will be effectively recovered by varying For example, by 
taking ~ 0.7 the Turk-Chieze and Lopez SSM [[J is reproduced. The factor is varied in 
the range 0.4 — 1.6 (the lower limit fs — 0.4 is in fact set by the Kamiokande result), which 
in fact corresponds to the variation of the central solar temperature within (—4 -f- +2) %. 

Let us first consider the vacuum oscillation solution in the case of two neutrino flavours: 
v x where u x is or v T (alternatively, it can be also a sterile state v s ). In this case 
the survival probability for solar u e 's with energy E depends on two parameters, the mixing 
angle 6 and the mass square difference 8m 2 : 

P(E) = 1 - sin 2 26 sin 2 (^^j (4) 

The sun-earth distance depends on time as L t = L[l — ecos(27rt/T)], where L = 1.5 ■ 10 n m, 
T = 1 yr, and e = 0.0167 is the ellipticity of the orbit. Then the time averaged detection 
rates expected in the radiochemical experiments read as 

R = j dEa(E) ^(P(E)4> k ) T Xk(E) (5) 

J k 

where cr(E) is the detection cross section, fc = fk<Po are the fluxes of the relevant components 
of the solar neutrinos (k = B, Be, etc.), Xk{E) are their energy spectra normalized to 1, and 
(. . .)t stands for the average over the whole time period T. For the Kamiokande detector 
we have 

R K = f dEX B (E) \(P(E)<p B ) T a Ue (E) + ((<p B ) T - (P(E)<p B ) T ) a Ux {E)\ (6) 

where a Vy (y = e, x) is the v y t~ scattering cross section (in the case of conversion into sterile 

neutrino a Vx = 0) and E t h = |(T e + JrjT e + 2m e )), where T e = 7.5 MeV is the recoil 
electron kinetic energy threshold. 

We accept the hypothesis that the solar neutrino luminosities are constant in time, and 
use the averaged data of the chlorine, gallium and Kamiokande experiments to perform the 
standard x 2 analysis for various cases. We find it instructive to separate the experimental 
and theoretical (SSM) uncertainties, and do not include the latter in \ 2 analysis once these 
are simulated by varying the factors We define 

X 2 = E(^T^) 2 (7) 

a V °a J 

where the index a = CI, Ga, K labels three different data and a a are the corresponding 
experimental errors. Once we describe the three experimental results by means of two 
paramemters (5m 2 and sin 2 29), the number of the degrees of freedom is 1. 

We have performed the \ 2 analysis for various values of Jb, using for the numerical 
minimization procedure the MINUIT package provided by the CERN program library. The 
corresponding best fit points of Xmin an d 'lc' areas containing the true parameter values 
with 68 % probability, once the solution is assumed, are shown in Fig. 1. We see that 
compared to the previous analysis [fLlp , the fits have become much better. E.g. in the case 



■ i 



/b = 1 we have Xmin = 1-6 (versus Xmin = 4-4 obtained in our previous paper JT4|), so that 
the just-so oscillation is allowed as a SNP solution at the 20 % confidence level (C.L.). [] 

By varying f B , the relevant range of 5m 2 remains rather stable, while sin 2 29 tends 
to be smaller with decreasing f B . The lowering (increasing) of /b implies the weakening 
(strengthening) of the oscillation. However, as a general tendency, by decreasing f B the fit 
becomes worse: e.g. for f B = 0.4 we have Xmin = 3-2, while it improves for f B > 1: e.g. for 
f B = 1.6 we have Xmin = 0-9' which is acceptable at 34 % C.L. 

Thus our conclusion is that for f B > 0.7 the relevant parameter range is limited by the 
values Sm 2 = (0.5 — 1) X 10~ 10 eV 2 and sin 2 29 = 0.7 — 1. For f B < 0.7 the parameter range 
is somewhat wider (see Fig. 1) but as we have seen in this case the quality of the fit is less 
good). The dependence on sin 2 29 is simpler - the just so scenario preferes the large (close 
to the maximal) mixing,^ whereas the dependence on Sm 2 is more sophisticated, due to 
the different oscillation length of solar neutrinos with different energies. To give some more 
numerical insight on this dependence, in the Figs. 2a,b we show x 2 as a function of Sm 2 for 
the cases sin 2 29 = 1 and 8/9, which correspond to the cases of the maximally symmetric 
(democratic) mixing in the 2 and 3 neutrino system (see below). 

Let us consider now the case of three neutrino mixing, when the neutrino mass eigenstates 
z/j with masses (i =, 1,2,3) are related to the weak eigenstates v a [a = e,/i, r) as v a = 
Uai^i, where U is the 3x3 unitary matrix of the neutrino mixing. The general expression for 
the u a — > up transition probability is: 



Pap{E) 



(8) 



The radiochemical detectors are sensitive only to v e state. On the other hand, the 
Kamiokande experiment does not distinguish between the contributions of and u T , pro- 
vided that the latter have no non-standard interactions with electrons {a v = o Vt in eq. @). 
Hence, the analysis of the solar neutrino data in the 3z/ case needs only the expression of the 
v e — > v e survival probability: 

P(E) = 1-4 \U el U e2 \ 2 sin 2 (^) - 4 \U el U e3 \ 2 sin 2 (^) - 4 \U e2 U e3 \ 2 sin 2 (^i) (9) 

4£/ 4i!/ 4£/ 

where Sij = m 2 — m 2 . This expression depends on four independent variables which we 
choose as U e 2,U e3 and S 2 i,6 3 i, whereas <5 32 = 5 3 i — 5 2 i and the unitarity constraints fix 
\U e i\ 2 = 1 — \U e2 \ 2 — \U e3 \ 2 . One can use e.g. Maiani-like parametrization by identifying 
U e i = cos 9 cos <fi, U e2 = sin 9 cos <p and U e3 = sin 0. 

The 2u oscillation is a limiting case which is achieved if the third neutrino is decoupled. 
E.g. in the limit U e3 = (sin0 = 0), in which case the dependence <53i is also lost, we have 
4:\U e iU e2 \ 2 = sin 2 29 and thus eq. (f|) reduces to the familiar expression @). On the other 
hand, when e.g. 5 32 = 0, we can identify £31 = 5 2 \ = 5m 2 and still 4|f/ e i| 2 (l — |[/ e i| 2 ) = sin 2 29. 

3 Since the theoretical uncertainties are not included, here and in the following the quantitative statements 
on the C.L. are somewhat informal, and reflect only the ideal situation with no uncertainties in SSM predicted 
values of the neutrino fluxes, exactly fixed by the choice of the factors fk- For the case of v x being a sterile 
neutrino we obtain Xmin = 5.5. 

4 According to ref. [ p^t , such a large mixing puts the just-so oscillation in contradiction with the observed 
pattern of the neutrino signal from SN 1987A. See, however, the recent paper Eq]. 
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In principle, the presence of four parameters versus only three experimental data makes 
the system overdetermined and thereby the standard meaning of \ 2 analysis is not directly 
applicable. However, even in this case the canonical \ 2 function fl7|) provides a numerical 
insight on the deviation of the theoretical estimates from the experimental data, within their 
errors. Its minimization with respect to all the four parameters can allow to know whether 
the 3u just-so scenario shows a better "fit" (i.e. smaller \ 2 ) lYi some relevant parameter 
region as compared to the 2v case. As smaller the value of x 2 is, as closer is the theoretical 
"just-so" prediction with experimental data. 

The result of our analysis is rather surprising. Performing numerically the minimization 
procedure in the four parameter space, we always find that the best fit point (for which \ 2 
achieves the minimum) is always located in the two parameter subspace which effectively 
corresponds to the 2v case, i.e. when e.g. U e s = or 5 32 = 0, for any feasible value of 
The values of the minimal Xmm coincides with that is obtained above in the 2v scheme. For 
example, in the case fs = 1 we obtain xLin = 1-6, which occurs for the parameter values 
\U e3 \ 2 = 0, \U e2 \ 2 = 0.40 (i.e. sin 2 29 = 0.96) and 5 21 = 0.6 • 10- 10 eV 2 . Of course, there 
are also trivially equivalent minima obtained by permutations U e i — * U e j, or by identyfying 
Oij = Sik- 

The feature seems rather obvious in the limit when say 631 is taken to be large, i.e. 
7713 ^> 10~ 5 eV. In this case the two last terms in (|9|) are in the regime of the averaged 
vacuum oscillation, and the v e survival probablity reads 

P{E) = 1 - - sin 2 2(f) - cos 4 sin 2 26 sin 2 (^|^) (10) 

Then for non-zero the data fit has to become worse, since the contribution of the energy 
dependent term becomes smaller which makes it more difficult to reconcile the Homestake 
and Kamiokande data. This feature, however, is not so obvious when 831 is also in the just-so 
(~ 10~ 10 eV 2 ) range. 

As we have seen, the best just-so description of solar neutrino data is achieved when 
in fact only 2 neutrino states participate the oscillation. Let us discuss briefly how this 
situation could be obtained in the context of various neutrino mass generation schemes. 

The neutrino mass range needed for the just-so oscillation can naturally emerge from the 



Planck scale effects |10|, [TT|. In the minimal standard model (or SU(5) GUT) the neutrinos 
are massless as far as renormalizable interactions are concerned. In other words, in the 
absence of the right-handed (RH) neutrino states the lepton number conservation arises as an 
accidental global symmetry of the theory due to the joint requirement of gauge invariance and 
renormalizability. However, the global symmetries need not be respected by nonperturbative 
quantum gravity effects. Thus, if the SM (or SU(5)) is a true theory up to the Planck scale, 
the neutrino masses can be induced only by the non-renormalizable operators cutoff by the 
Planck scale Mpf. 

-r^-{kH)C{l s K) (11) 

Mpi 

where k = {yi,ti) T {% = 1,2,3) are the left-handed (LH) lepton doublets, H is the standard 
Higgs doublet, and a^- are the order 1 constants. (In fact, this is equivalent to the seesaw 
mechanism with RH neutrino states having ~ Mpi Majorana masses and ~ 1 Yukawa 
constants; in particular, such a situation is automatically obtained in a supersymmetric 
grand unified model SU(6) fL9| ). As a result, the neutrinos get Majorana masses of the 



order of m = v 2 /M Pt = 3 • 1(T 6 eV, where v = 174 GeV is a VEV of the Higgs doublet. The 

parameter dm ~m thus naturally emerges in the range needed for the just-so oscillation. 

In general case, this operator induces mass terms for all three neutrinos u e , u^, v T . If 
the matrix a^- of the effective coupling constants is a general one, then we have the 3 — v 
oscillation case. However, it could occur that a^- is a rank-1 matrix, 

a = U T PU, P = diag(0, 0, 1) (12) 

where U is some unitary matrix. Then only one neutrino eigenstate gets mass of the order 
of m, while others stay massless. In this case the neutrino oscillation picture effectively 
reduces to the 2u case. For example, it was argued in that once the operator (|TT|) 
is induced by nonperturbative gravitational effects (virtual black holes or wormholes), the 
flavour blindness of gravity could imply the following pattern for the neutrino mass matrix 



U (13) 



in which case we effectively have 2v oscillation with 5m 2 = 9a 2 m and sin 2 26 = 8/9, in 
good agreement with the data fit in the just-so picture (see Fig. 2b). 

It is however also possible that one of the neutrino eigenstates acquires a mass larger 
than m from sources different from the Planck scale induced operators (|TT|). In particlular, 
the neutrino mass matrix may have a structure 

m v = mT + m", m 3>m (14) 
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where T is some rank-1 matrix. If T = P, where P is defined in eq. flT2|) , then the v T state 
acquires mass m and for m ~ 5 — 7 eV it can play a role of the cosmological hot dark matter 



component [21]. In this case for the SNP solution remains the just-so oscillation u e — > 
with maximal mixing sin 2 26 = 1, which still provides a good fit to the solar neutrino data 
(see Fig. 2a). 

Alternatively, one can imagine that the matrix T has a democratic structure in the space 
of the Up and v T weak eigenstates: 



T= Oil (15) 




In this case the — > v T oscillation emerges with sin 2 26 ^ T = 1 and <5m 2 T = 4m 2 , which for 
m ~ 5 • 10~ 3 eV could explain the the atmospheric neutrino deficit while the just-so 
oscillation is still effective in explaining the SNP. 

Concluding, we have seen that the just-so oscillation scenario provides a very good fit to 
the recently updated experimental data on the solar neutrinos, and its present status is as 
reasonable as that of the widely popular MSW scenario. We have shown that best data fit 
is achieved within the two neutrino system, and with incorporating the third neutrino the 
fit can only get worse. This points that the third neutrino state can be involved into the 
games of solving other present neutrino puzzles as are the atmospheric neutrino problem or 
the problem of the cosmological hot dark matter. 
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Figure Captions 



Fig. 1. The long-dashed, dott-dashed, solid, dashed and dotted contours mark the 68 % 
C.L. regions respectively for f B —0A , 0.7, 1, 1.3, 1.6, The corresponding best fit points are 
shown by diamonds, with Xmin being respectively 3.2, 2.5, 1.6, 1.5, 0.9. 

Fig. 2. x 2 as a function of Sm 2 for a fixed mixing angle: sin 2 29 = 1 (a) and sin 2 29 = 8/9 
(b). The dot-dashed, solid and dashed curves correspond to the cases fs = 0.7, 1 and 1.3, 
respectively. 

Fig. 3. x 2 minimized with respect to parameters sin 2 29 and S21, as a function of the 
'13' mixing sin 2 20, for the case of large S 31 (m 3 3>m). The dot-dashed, solid and dashed 
curves correspond to the cases fs = 0.7, 1 and 1.3, respectively. 



